Electrostrictive counter-force on fluid microdroplet in short laser pulse 
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Compiled March 1, 2013 
When a micrometer-sized fluid droplet is illuminated by a laser pulse, there is a fundamental distinction 
between two cases. If the pulse is short in comparison with the transit time for sound across the droplet, 
the disruptive optical Abraham-Minkowski radiation force is countered by electrostriction and the net stress 
is compressive. In contrast, if the pulse is long on this scale, electrostriction is cancelled by elastic pressure 
and the surviving term of the electromagnetic force, the Abraham-Minkowski force, is disruptive and deforms 
the droplet. Ultrashort laser pulses are routinely used in modern experiments, and impressive progress has 
moreover been made on laser manipulation of liquid surfaces in recent times, making a theory for combining 
the two pertinent. We analyze the electrostrictive contribution analytically and numerically for a spherical 
droplet. © 2013 Optical Society of America 



Consider a laser pulse impinging on a microdroplet 
of a homogeneous fluid, whose radius is a and of den- 
sity p, situated at the origin. Assume the pulse to be 
short (defined below), but long enough to be treated as 
a plane monochromatic wave. The refractive index of the 
droplet is n, its permeability /.t = 1, and we assume it is 
surrounded by vacuum for simplicity. The fluid is then 
acted upon by two different kinds of optical forces; a 
volume force in the interior due to electrostriction (ES) 
acting towards regions of higher electromagnetic energy 
density, and a surface force acting on the surface r = a 
due to the difference in permittivity between fluid and 
vacuum. 

When the pulse length is sufficiently long, the ES vol- 
ume force is compensated by mechanical fluid pressure, 
and the fluid motion becomes entirely dictated by the 
surface force and fluid mechanical equations of motion. 
This was the case in the classical experiment by Zhang 
and Chang [1], and to our knowledge only this case has 
been considered in the theoretical literature [2-5]. The 
time it takes for the electrostrictive compressive force to 
be compensated depends on the compressibility of the 
fluid; it is approximately the time it takes for a sound 
wave to traverse the droplet, Tc = 2a/u where u is the 
speed of sound. 

The effect has become accentuated in recent years be- 
cause of the frequent use of short laser pulses. Moreover, 
in some modern applications one works with very small 
surface tensions, typically in a two-fluid system such as 
used by Delville's group, c.f. [6] and references therein, 
in which surface tension is reduced to a millionth that 
of an air-water surface. The surface is thus sensitive to 
even small stresses. 

The electromagnetic force density in a fluid is [7] , 
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X> = eoe£, 3 ~ fJ-o'H] thus e is nondimensional, and 
we assume it to be real and positive. We neglect an ad- 
ditional Abraham term (whose existence is subject to 
debate), which is much smaller than the other in all 
typical parameter ranges [7]. Calligraphic typeface de- 
notes real held quantities, so £ = Re{i?e*'^*}, etc. Note 
(f^) = ^|i?P and so on. 

Consider flrst the second term in Eq. (1), the elec- 
trostriction term, called f^^. We draw on the Lorentz- 
Lorenz relation to evaluate the derivative of the permit- 
tivity, whereby 
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One sees that the effect from this force on the droplet 
is compressive, as it points towards higher optical in- 
tensity. The striction force, when taken alone, acceler- 
ates the fluid inward until an elastic counter pressure is 
established, on a time scale Tc- We recently considered 
pressure waves set in motion by a laser beam in a homo- 
geneous fluid [8]. 

Secondly, there is the flrst term in Eq. (1) which may 
be called the Abraham-Minkowski (AM) term f^^^^ as 
it is common for the Abraham and Minkowski energy- 
momentum tensors [7] . This term is the only one to sur- 
vive after the elastic response time as discussed above. 
Thus it is sufficient to describe the laser-induced sur- 
face deformations observed in quasistatic or long pulse 
experiments, f.ex. Refs. [1,6,9]. 

By integrating the force density (2) across the bound- 
ary we obtain the ES surface pressure P^^ , 
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where (• ■ • } denotes average over an optical period. Ac- 
cording to our conventions the constitutive relations are 



where n is the outward normal and cr^^ the scalar pres- 
sure, a negative quantity. In the following we shall as- 
sume circularly polarized plane wave illumination, draw- 
ing on the Mie scattering formulation of Barton et al. 
[10]; the circularly polarized fields are given in full in [5]. 



According to Eq. (2) the ES force also has a volume 
contribution from the interior of the sphere proportional 
to V(S^. One could calculate this force in detail, yet for 
our purposes it is physically more instructive to consider 
some overall properties of the ES volume force which 
allows simple comparison to the AM surface force. 

First, note that the net ES force on the whole sphere 
is zero. This is easily seen: the force (2) is a gradient 
and when integrated over any volume V can be written 
as an integral over surface area. For any control volume 
outside the sphere, de/dp is zero, hence zero total force. 

Secondly, a physically instructive quantity to calculate 
is the full ES force (surface and volume contributions) 
acting on the front and back hemispheres. Choosing a 
control volume bisecting the sphere at ^ = and closed 
outside the sphere for z < 0, the force may be written 
as an integral over the circular section only. 
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r being the radius in cylindrical coordinates. < (>) de 
notes front (back) as seen by light propagating from z 
— cxD. Exactly the same argument gives Ff^ = 
which accords with zero total force. 

Consider next the AM force f^^, which according to 
Eq. (1) acts in the dielectric boundary layer only, where 
Ve ^ 0. The corresponding surface pressure P^'^ is 
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Here £t = £e(^ + £(t>4> is the field component parallel 
to the surface, and £r the component normal to it. The 
expression (5) is positive, in accordance with the general 
property of the dielectric gradient force that it always 
acts in the direction of the optically thinner medium. 

For comparison we calculate the AM force on the front 
(<) and back (>) hemispheres. 
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and similarly for the back upon replacing the integral 
with /„ . Unlike for ES, the two arc of unequal magni- 
tude and the net force is a push in the direction of light 
propagation; the sphere acts as a lens concentrating the 
electromagnetic energy near the rear surface. 

The field equations for circularly polarized light are 
straightforward to work out in the formalism of [10], 
and are quoted, e.g., in Appendix C of [5]. We let 
a = 27ra/A_ = nuja/c be the number of circumfer- 
ences per wavelength at frequency lo. With field com- 
ponents Er, Eg, E^ written out, the front and back 
hcnnsphcre AM forces can be calculated numerically, 
as shown in figure 1 where we plot dimensionless force 

gAM ^ ^AM/(^^2g^^2)^ ^j^g^g ^2 -g ^j^g ^.^^g ^^i^^ ^f 

the (real) electric field. 

We are able to calculate the hemisphere total ES 
forces analytically, since at 9 = ^ Lcgendre polynomi- 
als simphfy through P^(0) = ~Sni, -F',1'(0) = Sno and 



AM force reorfiolf . 
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Fig. 1. Force components on front and rear hemispheres. 
Q is repulsive if positive for rear half or negative for front 
half. 
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sin a; — cosx. After some straightforward 



calculation we find F^% = ±TTa'^eoE^Q^^ with 
gES =2(?i2 _ i)(n2 + 2){|cin4/i(na) + hina)] 
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+ x^ (2a; cos a; — sin x) sin x'] 
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Here /i.2,3 are the radial integrals L dx over ijJi^k)/ k^ , 
iPI{x)/k and 'il)'i{x)/x, respectively, resulting from in- 
sertion of E'^ into Eq. (4) . Ci is the cosine integral. 

In figure 1 wc compare the hemisphere forces on the 
front and back hemispheres due to electrostriction and 
electromagnetic forces. In all figures n = 1.33. The to- 
tal AM force, dictating the centre of mass motion of the 
droplet, is positive (pushed by the laser), whereas the to- 
tal surface force is much smaller and of of opposite sign, 
and the total surface force slightly is compressive near 
the focal points at the droplct's rear. ES volume forces 
balance the corresponding surface forces in sum since 
also the droplct's interior liquid is attracted to the areas 
of large field strength near, but not at, the droplct's rear 
surface. According to Eq. 2 {£'') acts as a potential for 
the electrostrictive force density, and is shown in Fig (2) 
for four values of a. The higher a, the more electrostric- 
tive compression is localised in small areas near the rear 
droplet boundary. 

In Fig. 3 we show how optical surface forces are distrib- 
uted over polar angles 6*, plotting (j^^{9), cr^^(cr) and 
the sum of these, (T*°* . The net surface force is compres- 
sive, which serves to counteract the disruptive AM force: 
The AM surface force is nearly cancelled by the ES com- 
pressive force, leaving a much smaller cr*°* which is com- 
pressive on average but can take small repulsive values in 
certain areas. The cancellation is exact for n — > 1, while 
the net surface force density becomes more strongly com- 
pressive for increasing n. Note correspondence with first 
two panels of Fig. 2. 




Fig. 2. Elcctrostiction potential {S^)/E^ inside droplet 
for four values of a. Laser beam incident from left. 
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Fig. 3. Surface force densities a = a/{e^E'^) for AM 
and ES as functions of polar angle 6. Polar plot of r = 
&{9) + 20 for illustration (20 is an arbitrary number for 
visualisation). 



Wc have calculated the effect of clectrostriction for 
the optical force on a liquid droplet. The resulting be- 
haviour depends crucially on the laser pulse duration 
compared to the sonic transit time across the dimensions 
of the sphere. For example, in water at room tempera- 
ture, u K, 1500 m/s, so that if a = 50 /im, the pulse is 
short in this sense when t < 70 ns, which is well achiev- 
able in practice (current ultrashort pulses have femtosec- 
ond duration). For the purpose of optical pulling [11] of 



microdroplets, pulses of a few ns at intervals of about one 
//s were found to be appropriate [5]. Laser pulses used 
in droplet manipulation experiments have convention- 
ally been longer than this, allowing the elastic counter 
pressure time to build, cancelling the contribution from 
clectrostriction. Only the Abraham-Minkowski force re- 
mains, and the force is repulsive as Fig. 1 shows. This 
is sufficient to explain experiments conducted to date, 
e.g. [1,6]. At times much earlier than t <C 2a/u, how- 
ever, the liquid respons both to the electrostrictive and 
AM forces as there is no appreciable counter force met 
by elastic pressures (there is in this sense an analogy 
with the classic non-equilibrium pressure experiment of 
Goetz and Zahn from around 1960 [12]). Inclusion of 
clectrostriction is now crucial to describing the behaviour 
of the droplet. 

In order to describe droplet hydrodynamics at short 
times, we have shown that a fully compressible theory is 
required, whereas previous theoretical treatments were 
incompressible [2-5]. What the effect on the droplet 's 
motion from clectrostriction will be is still and open 
question of great interest. No experiments exist to our 
knowledge where this has been probed. It is clear that 
the initial deformation is to compress the droplet. If vis- 
cosity is low, overcompensation by the elastic pressure 
build-up would be expected to give rise after a while to 
surface deformations similar to those for a longer pulse, 
whereas with higher viscosity we expect surface oscilla- 
tions to be significantly smaller due to the conteraction 
of clectrostriction. It seems to us that the development 
of a fully compressible theoretical framework for this ge- 
ometry is a natural goal for the near future. 
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